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A technique to accomplish fu //-spectrum arraying where all the telemetry power 
is put into the subcarrier sidebands (suppressed carrier,) is described. The matched 
Biter needed in each antenna prior to cross-correlation for deriving the coherence 
delay and phase offsets is an open-loop version of the telemetry phase-lock loop 
provided in the Advanced Digital Receiver . In analogy with a Costas-loop teleme- 
try receiver , a “squaring loss ,} is derived and a signal-to-noise ratio for the cross- 
correlation loop phase is presented. 


I. Introduction 

Normally, as a spacecraft travels farther from Earth 
and the telemetry signal-to-noise ratio (SNR) gets poorer, 
two system parameter trade-offs come into play. First, 
the telemetry modulation index is usually increased so 
more transmitter power is moved from the carrier to the 
telemetry signal, thereby improving telemetry SNR. This, 
of course, may result in a carrier signal that is significantly 
harder to acquire and track. The limit for this trade-off 
is full modulation where no carrier power is present. In 
this case, the carrier signal frequency must be acquired 
and tracked using a less-than-optimal Costas phase-lock- 
loop technique. The capability to Costas-loop track is not 
presently available in the DSN, but is planned as part of 
the new Block V receivers. 

The second trade-off that comes into play is the rate 
at which telemetry data are transmitted back to Earth; 
this rate can be reduced, resulting in an improved SNR 
per telemetry bit. This has the unfortunate consequence 
of also reducing the total amount of data that can be re- 
turned during the critical encounter phase of a mission. 

A technique often applied within the DSN to overcome 
these constraints (short of building larger antennas, even 


lower noise receivers, or employing more advanced data 
encoding/decoding methods) is antenna arraying to en- 
hance the SNR of the received telemetry [1]. Perhaps the 
most recent major applications of antenna arraying oc- 
curred during the Voyager 2 encounter with Uranus and 
the Voyager 2 encounter with Neptune. In each case, ar- 
raying provided adequate telemetry SNR employing data 
rates higher than would have been possible otherwise. 

At least four different techniques, which depend on the 
details of the observing circumstance, have been used for 
arraying or combining the signals from several antennas. 
These include symbol-stream combining, baseband array- 
ing, carrier arraying (or more correctly, carrier aiding), 
and full spectrum combining. In the various discussions 
of these techniques, the terms arraying and combining are 
usually used interchangeably. 

The first technique, symbol-stream combining, works 
well when each ground-based antenna/receiver is capable 
of locking on the spacecraft telemetry stream and demodu- 
lating it down to the soft symbol stream (the raw telemetry 
data before the decision of whether a given bit is 1 or 0 
is made). Each stream must then be delayed to line up 
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symbols, followed by weighting and combining. One major 
advantage of this technique is that the data stream from 
each antenna arrives at only the telemetry rate (typically 
on the order of 100 kbaud), facilitating its transmission to 
the combining location, either in real time or through a 
recording medium. The International Cometary Explorer 
(ICE) mission successfully used symbol-stream combining 
in a primary mode to gain approximately 2 dB in SNR 

[2] . A short time later, this combining technique was used 
on Voyager 2 to gain slightly less than 3 dB improvement 

[3] . More recently, symbol-stream combining was used as 
a backup for baseband arraying during the Voyager 2 en- 
counter with Neptune [4]. 

The second technique, baseband arraying, works when 
the telemetry SNR is high enough to permit locking on 
and tracking of the carrier signal at each antenna, but 
too low to reliably maintain subcarrier lock and symbol 
sync. In this case, the baseband signal obtained after 
carrier demodulation is delayed, weighted, and added to 
a corresponding baseband signal from the other antenna. 
Then subcarrier demodulation and symbol sync are ac- 
complished on the combined stream. The baseband band- 
width is on the order of 3 MHz and therefore somewhat 
more difficult to transmit or record for combining. During 
the Voyager 2 encounter with Uranus, this technique was 
used effectively with the Parkes Radio Astronomy antenna 
in Australia [5]. 

The third technique, carrier aiding, is useful when one 
antenna in the array can acquire and maintain carrier lock, 
but the other antennas cannot, either because they are 
smaller or have higher noise receivers. If the antennas to 
be arrayed are in proximity, the lock signal from the more 
sensitive antenna can be sent to the other antennas in real 
time and used to aid their locking on the telemetry carrier. 
The resulting baseband signals are then treated as they 
are in baseband arraying. This technique has been used 
quite recently to successfully array a 70-m antenna and a 
34-m antenna at Goldstone while looking at the Pioneer 11 
spacecraft. 1 

Finally, full-spectrum combining of an open-loop space- 
craft telemetry signal can be used when the carrier is too 
weak to track, or when it is not convenient to track, at 
any single antenna. One form of this technique was used 
to coherently add the signal from all the antennas of the 
Very Large Array (VLA) during the Voyager 2 Neptune 


1 T. Peng, “Carrier Array Demonstration at Goldstone with Pio- 
neer 11,” JPL Interoffice Memorandum 3303-90-62 (internal docu- 
ment), Jet Propulsion Laboratory, Pasadena, California, May 11, 
1990. 


encounter [6]. Following is a description of a significantly 
improved version of full-spectrum combining, with com- 
ments on how it differs from the original. 

By way of introduction, note that the telemetry signal 
transmitted from a spacecraft can be represented by the 
equation [7] 

S[t] = \/2P{cos A sin [Q c f] + d[t ] sin AEs*[<]} 
where each term in the summation (odd k from 1 to oo) is 
s k [t] = (2/tt)(1 /k){g uk sin [( Q c + kQ 9C )t] 

- gi k sin[(Q c - fcfUc)*]} 

corresponding to the subcarrier sidebands. The amplitude 
corresponds to a signal total power of P. Figure 1 is a 
graph of this signal. 

In these expressions, 
t — time 

Q c = the transmitter carrier frequency 
Q sc = the telemetry subcarrier frequency 
A = the modulation index of the subcarrier 
k — the subcarrier harmonic number 
d[£] — the telemetry modulation 
g uk — the upper sideband gain at frequency k 
g lk — the lower sideband gain at frequency —k 
P = the signal total average power 

For simplicity, assume the modulation index is 90 deg, 
so that only the suppressed carrier case is treated (this will 
likely be true when full-spectrum combining is used to put 
all possible power into the telemetry). Also, assume the 
gain factors are all unity so the summation over k can be 
treated as a square wave (the efFects of this on the end 
result are minor). The signal then becomes 

S[t] = y/2 Pd[t] cos [Q c t)S[Q sc t] (1) 

where the outlined S represents a square wave with zero 
crossings coincident with a sine wave of the same argu- 
ment. 
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The spacecraft signal in Eq. (1) is received by several 
antennas, located at various delays, r, from the spacecraft. 
These delays include geometric as well as nondispersive 
media and instrument components (as well as a large com- 
ponent corresponding to the distance from the spacecraft 
to the center of the Earth and common to all antennas, 
which will be ignored because it differences out in the fi- 
nal results). In addition, each antenna’s receiver down- 
converts the signal to a carrier intermediate frequency, 
Qif = (ft c — using a local oscillator frequency, 

The combined effect on the signal received at the mth an- 
tenna is 

5 m [t] = V2Pd[t - r m ] cos [fy/i - Cl c r m + 0 om ] 

X S[a, c (t-T m )] + N m (2) 

where 9 om is an unknown, slowly varying (its change in a 
time equal to the difference in the delays between the var- 
ious antennas is negligible) residual phase due to various 
dispersive media and instrument effects. N m is the noise 
of variance <r om in the bandpass containing the signal. 

The goal of the combining technique is to coherently 
add this signal from several antennas to obtain an im- 
proved SNR for telemetry extraction. To achieve this goal, 
the received signals must be delayed and phase shifted to 
correct for the relative effects between the various antennas 
that are represented in Eq. (2) before they can be added. 
In the following analysis, assume the signal (plus noise) 
modeled in Eq. (2) is digitally sampled at the Nyquist rate 
and processed with accuracy sufficient to ignore quantiza- 
tion and round-off effects. 


II. Delay and Phase Shift 

Following the steps outlined in Fig. 2, begin by delaying 
the signal received at the mth antenna by t' mi based on a 
best a priori model, r f my of the antenna/spacecraft system 
(the caret on this delay signifies a quantized version of the 
model delay; the prime indicates a model delay). From 
this, 

= \/2 Pd[t - A m ] cos + f' m ) - Cl c T m + 6 om ] 

x S[Q JC (< — A m )] (3) 

where the quantity A m = r m — r f m is the residual between 
the actual delay and the value used in data processing. 


The delay operation is followed by a phase shift that 
involves first the generation of a quadrature version of 
Eq. (3) using a hybrid and then a complex multiplication 
by an expression of the form 

Epm[t] = exp [-jOpm] (4) 

The final signal becomes 

Sm[<] = V2Pd[t - A m ]exp[?'{n i/ (< + f^) - fi c r m + 0 m }] 

xS[fi JC (<- A m )] (5) 

where 9 m — 9 om — 9 pm is the resulting phase offset of the 
signal. Note that the signal is now a complex quantity, 
having both in-phase and quadrature-phase components. 

By adjusting T* m and 9 pm , it is possible to bring the 
signals from each antenna into coherence for combining. 
An error signal will be generated by comparing the sig- 
nal in Eq. (5) between pairs of antennas and then using 
this signal to control the delay and phase adjustment of 
the telemetry streams from these antennas relative to each 
other. 


III. Matched Filter 

The first step toward obtaining an error signal is to 
perform a matched filtering operation on the signal. When 
implementing full spectrum arraying with the VLA during 
the Voyager Neptune encounter, this “matched” filtering 
consisted simply of narrowing the VLA receiver bandpass 
from its nominal 50 MIIz down to about 8 MHz with base- 
band filters. While hardly optimal, such a procedure is ac- 
ceptable if the SNR at each antenna is reasonable, as was 
the case for Voyager. However, in general one needs a more 
optimal approach. A true matched filtering (providing a 
better power SNR by as much as a factor of 8 MIIz/43 kHz, 
or ~ 180 in this case) is properly accomplished through the 
implementation of two procedures: (1) carrier demodula- 
tion followed by (2) sub carrier demodulation. 

A. Carrier Demodulation 

Carrier demodulation is accomplished by coherently de- 
tecting the signal at the carrier intermediate frequency 
with a locally generated signal of the form 

E cm [<] = exp [- j + f^) - Sl' c T' m + 0 em }] (6) 
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where the primed quantities are the best guesses for their 
unprimed counterparts, and 6 cm is an unknown arbitrary 
phase. The result of this process is a signal with real 
and imaginary parts representing the in-phase (I) and 
quadrature-phase ( Q ) components of the demodulated car- 
rier. After low-pass filtering, 

< I cm > = (1/2)72 Pd[t - A m ] cos [$ cm ]S[ft ,c{t ~ Am)] 

(7) 

< Qcm > = {l/2)V2Pd[t - Am] sin [$cm]S[ft JC (< - Am)] 

( 8 ) 


where 

<$ cm = (fi c - Q' c )t - 4- 4 - & cm (9) 

is the residual carrier demodulation phase and, with good 
modeling, is only a slowly varying function of time. In 
Eq. (9), the quantity A' m = r f m — is the model residual 
delay. 

B. Subcarrier Demodulation 

Subcarrier demodulation is accomplished by coherently 
detecting the signals in Eqs. (7) and (8) at the subcarrier 
frequency using 

E, m [t] = jEH{fUf - Am) + <?,m}] (10) 

where ft' c is the best guess of the subcarrier-oscillator fre- 
quency, 9, m is an unknown arbitrary phase, and the out- 
lined E represents a complex square wave with components 
c + j S. Again, the results of this process are signals with 
real and imaginary parts representing the in-phase and 
quadrature-phase components of this subcarrier demodu- 
lation. After filtering, there are four components: 

< IcmI.m > = (l/2)%/2 Pd[t - Am] 

X COS [<f>cm]C + [$ $m ] + A him (11) 

< QcmI.m > = (1/2 )V2Pd[t - Am] 

x sin [<*cm]C + [<I> jm ] + N QIm (12) 


< IcmQ am > = (1/2)72 Pd[t - Am] 

X COS [$ cm ]S + [$ J m] + NjQm (13) 

< QcmQ.m > = (1/2)72 Pd[t - Am] 

X Sin [<J> cm ]S + [$,„,] + NQQm (l 4 ) 

where 

$$m = (&SC - fyc)* “ 

is the residual (slowly varying) subcarrier demodulation 
phase. Specifically, the real component of Eq. (10) when 
mixed with Eqs. (7) and (8) and then filtered gives 
Eqs. (11) and (12); the imaginary component of Eq. (10) 
when mixed with Eqs. (7) and (8) and then filtered gives 
Eqs. (13) and (14). The outlined C + and S+ quanti- 
ties represent the convolution of appropriate square waves 
(sine or cosine) as a function of their offset in phase, these 
waves have been lowpass filtered. They are, in fact, trian- 
gle waves as a function of this phase argument, with peak 
amplitude 1. The four independent noises, AT/jm> 

NiQ m , and N QQm , have been included explicitly at this 
point, each having a mean square variance of 

O-m = (l/4)^m/ n a ( 16 ) 

In this equation, recall that <r m is the noise deviation 
per data point (noise coherence interval) as represented 
by N m in Eq. (2), while the 1/4 factor is due to averaging 
over the IF cycles and n a is the number of data points 
over which the lowpass filters average. It is desirable that 
the value of n a correspond to a length of time equivalent 
to a telemetry bit interval (the largest it can be when the 
telemetry, represented by the function d(t) } is unknown). 
The goal of the delay and phase modeling is to minimize 
the changes in residual delay and phase so this integration 
can be as long as possible before the cross-correlation is 
performed. Note that to achieve an n a corresponding to 
a full telemetry bit interval, the symbol-sync on each of 
the data streams must be obtained independently. This 
is possible in the face of low SNR, since one can integrate 
over many telemetry bits to acquire sync. 

IV. Cross-Correlation 

The final step in the detection process requires the 
“squaring” of the signal to eliminate the telemetry data, 
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d[t] y which has an amplitude of ±1 but an unknown code. 
This is accomplished by cross-multiplying, or correlating, 
the signals from pairs of antennas. Unfortunately, this is a 
noisy process because of the multiplying of noise by noise, 
and noise by signal that comes about when generating the 
desired signal-by-signal product. If the SNRi of the incom- 
ing data is low (<1), then the SNR 0 of the outgoing signal 
will be a function of the square of SNRi , which deterio- 
rates rapidly with decreasing SNRi. The loss of SNR due 
to this effect in a Costas loop is usually called “squaring 
loss.” Here it more appropriately should be called “cross- 
correlation loss.” 

By taking some of the possible products (traceable from 
their names) of Eqs. (11) through (14) pair-wise for two 
antennas (the mth and the nth), for the signal-by-signal 
part, 

< IlmQIn > = ( P/2)d + [|A m - A„|] COS [$ cm ] 

x sin [$„] C + [S, m ] C + [$,„] (17) 

< QI m II n > = ( P/2)d+ [|A m - A„|] sin[$ em ] 

X COS [<£ C n] C + [$, m ] C + [$„,] (18) 

< IQmQQr, > = (P/2)d + [|A m - A„ |] cos[$ cm ] 

x sin [$ c „] S + [$ jm ] S + [$,„] (19) 

< QQmIQn > = ( P/2)d + [|A m - A„ |] sin [<S> cm ] 

x cos [$ cn ] S + S + [$ s „] (20) 

where 

d+ [|A m — A n |] = < d[t — A m ]d[t — A n ] > (21) 

is the time average of a convolution of the telemetry data 
as a function of the delay residuals. This function has a 
maximum value of 1 at the origin, and drops linearly to 
zero at an absolute delay difference of one telemetry bit 
length, assuming d[f] takes on values of ±1, and successive 
data bits are uncorrelated. 

If Eq. (17) is subtracted from Eq. (18), and Eq. (19) is 
subtracted from Eq. (20), 


Qcc = ( P W d+ [|Am - An |] Sin [<*> cm - <D C n] 

X C + [$ 3m ] C + [$,„] 

Qbs = (e/2)rf + [|A m - A„|] sin [$ cm - $ cn ] 

X S + [*, ro ] S + [$,n] 

When these two expressions are added, 

Qc = ( p /W+ [|Am - AJ] sin [$ cm - <f> cn ] 

X C + [<f> sm - $ jn ] (22) 

In this expression, C + [$, m - $ 5n ] is the time average of 
the sum of C + [^ sm ]C _|1 [$,„] and S + [$ 5 m ]S + [$ sn ] and is 
simply a function of the difference of and $ sn . In fact, 
this function is a slightly smoothed version of the triangle 
function represented by C + [$] (hence, the caret on C + ) 
with a maximum value of 2/3 at $ = 0. 

By taking all other possible products of Eqs. (17) 
through (20) for two antennas and combining these in a 
fashion analogous to that used to derive Eq. (22), 

h = (P/2)d + [I Am - A„0 cos [$ cm - $ cn ] 

X C + - $,„] (23) 

Qa - ( P/2)d + [|A m - A„ |] sin [$ cm - <f> cn ] 

x S + [<f 3m - $ , n ] (24) 

7 s = ( 7 7 2 ) <7+ (I Am - A„|] cos [$ cm - $ cn ] 

x S + [$ 3m - $ sn ] (25) 

The total cross-correlation noise on each of the compo- 
nents given in Eqs. (22) through (25) is due to the sum 
of the products of the various signals and noises explicitly 
shown in Eqs. (11) through (14). It is straightforward, but 
tedious, to show that these cross-correlation noise compo- 
nents are all independent of each other, have zero mean 
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value, and all have the same maximum mean square vari- 
ance of 

a 2 — 4cr 2 (P/6 + cr 2 )/n t (26) 

In Eq. (26), P is the power in the original telemetry signal, 
a 2 is the mean square noise variance per telemetry bit for 
the two uncorrelated signals making up this cross product 
from Eq. (16), 

ff2 = + < r n )/2 ( 27 ) 

and n t is the number of telemetry bits over which the 
cross- cor relation sum is averaged. 

V. Residual Delay and Phase Estimation 

Following the practice used in very long baseline inter- 
ferometry (VLBI) cross-correlation, each of these signals 
can be determined for several delay lags (i.e., for several 
values of |A m — A n | around the best a priori estimate). 
Then, based on the signal amplitudes in each lag, the cor- 
rect delay difference can be determined and used to adjust 
the delay of one of the data streams relative to the other 
to maintain data stream delay coherency. In actual fact, 
when the telemetry bits are many microseconds in length, 
the delay difference can be estimated a priori to a small 
fraction of this bit length, eliminating the need for multilag 
correlation. 

The phase differences in Eqs. (22) through (25) can be 
written explicitly as 

m - = ^(A^ - An) - Q c (A m - A n ) 

+ $ om @on @pm 

“h @pn ~~ @cm + @cn ( 28 ) 

and 

= ^c(^m — A' n ) — 0, c (A m — A n ) 

— C + (29) 

where A* = r k - r' k and A' k = r k - f' k are the residual 
delay and the model residual delay for the A:th antenna. 


Assuming (1) the frequencies can be specified to within 
a few hertz, (2) the rate-of- change of the delay differen- 
tials is no greater than a few picoseconds per hour (quite 
easily obtained, based on VLBI experience), and (3) the 
instrument phases are indeed slowy varying (a few tens of 
millihertz), it is possible to estimate these phase differences 
and feed them back to the appropriate phase mixers in one 
of the antenna data streams to maintain data stream phase 
coherency. Note that for the carrier phase in Eq. (28), this 
correction can be accomplished in any one of three phase 
mixers, while for the subcarrier phase in Eq. (29), it must 
be performed in the subcarrier demodulation mixer. 

In actual practice, the subcarrier phase residual in 
Eq. (29) can be kept small. This is because the delay resid- 
uals and the antenna clock offsets can usually be main- 
tained at less than 25 nsec, which, for a subcarrier fre- 
quency on the order of 1 MIIz, will result in no more than 
10 deg of phase error. Therefore, the only residual that 
need be estimated is that for the carrier phase. 

VI. Signal-to-Noise Ratio 

Estimating the noise in measuring the cross-correlation 
phase, as represented in Eq. (28), is completely analo- 
gous to determining the same quantity for a Costas phase- 
lock loop when tracking telemetry. As noted above, the 
Costas-loop algorithm involves the multiplying of quadra- 
ture signal components from a single antenna, while cross- 
correlation involves multiplying analogous signal compo- 
nents, but crosswise from a pair of antennas. 

Using the quantities given in Eqs. (11) and (12) from 
the same (e.g., mth) antenna, the Costas-loop equivalent 
of Eq. (22) can be derived: 

Q c = ( P/4)d + [0 ] sin [2<K cm ] C + [$. m ] (30) 

In an analogy with Eq. (26), the mean square variance for 
this Costas component is 

<^c = (r m ( P / 2 + a m)/ n t ( 31 ) 

From these two equations, the signal- to-noise ratio for 
measuring the carrier phase can be derived: 

SNR C = p/{ 1 + 1/(2 R d )} = pS L c (32) 

where p = 2P n a n t /<r 2 0 is the equivalent SNR in the final 
loop bandwidth, Rd — P n a /cr 2 is the data (symbol) SNR, 
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and Sic is the squaring loss as a function of Rj due to 
multiplication of signal by noise, and noise by noise, in 
the Costas- loop process: 

S u [Rd] = l/{l + l/(2^)} (33) 

Following an analogous procedure using Eqs. (22) and (26), 
the SNR for measuring the cross-correlation phase can be 
derived: 

SNR X = p (1/3)[1/{1 + l/(27?d/3)}] = pS Lx (34) 

The squaring loss here is the same as that for the Costas 
loop, except for scaling: 

S Lx [Rd} = (l/3)5xc[(l/3)i?d] (35) 

The larger loss in the cross-correlation case is due to the 

fact that the function C in Eq. (22) has a peak value of 
2/3, while the function C + in Eq. (30) has a peak value 
of 1. The only way to overcome this loss is by making 
nt, the number of symbols over which integration takes 
place, larger. This is analogous to making the Costas loop 
bandwidth narrower. 

Figure 3 presents plots of these two squaring losses as 
a function of Rj. It is evident that to obtain an SNR for 
the cross-correlation phase equal to that for the Costa s- 
loop phase, integration must be at least eight times longer 
when Rd is less than 1. 


Because of instability in the frequency of the spacecraft 
transmitter, increasing the integration time (or equiva- 
lently, narrowing the loop bandwidth) of the Costas-loop 
case beyond about 1 sec is not possible without serious 
SNR loss. But for the cross-correlation case, this insta- 
bility largely “common-modes out” from one antenna to 
the other, and therefore integration periods up to 100 sec 
are easily obtained with only slight SNR degradation. The 
limit on this integration time is set by the differential tro- 
pospheric effects between the two antennas. 

VII. Conclusions 

A technique has been outlined for performing full- 
spectrum combining of telemetry signals. It is possible 
to implement a device to accomplish this type of array- 
ing by using the DSN Advanced Receiver (ARX). ARX 
components would serve as the matched filters while new 
hardware would have to be added to provide the incoming 
delay and phase rotation. The real practicality for the use 
of full-spectrum combining in the DSN ultimately rests on 
the difficulty and expense of implementing these matched 
filters. It is expected that continued advancement in VLSI 
technology will easily provide inexpensive solutions in the 
near future. 

An alternative approach for a demonstration of the ba- 
sic ideas would be to record telemetry data with the DSN 
wide channel bandwidth (WCB) VLBI system and use the 
Block II VLBI correlator to accomplish the combining. 
Only a slight modification in hardware together with some 
software upgrades would be needed to test the technique. 
The one-bit sampling of this system, however, would limit 
its value as a combiner in any real application. 
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Fig. 1. Spacecraft telemetry signal. 
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Fig. 2. Delay and phase shift. 
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Fig. 3. Squaring loss. 
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